Friction and diffusion of matter-wave bright solitons 
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We consider the motion of a matter-wave bright soliton under the influence of a cloud of thermal 
particles. In the ideal one-dimensional system, the scattering process of the quasiparticles with the 
soliton is reflectionless, however, the quasiparticles acquire a phase shift. In the realistic system of 
a Bose-Einstein condensate confined in a tight waveguide trap, the transverse degrees of freedom 
generate an extra but small nonlinearity in the system which gives rise to finite reflection and leads 
to dissipative motion of the soliton. We calculate the velocity and temperature-dependent frictional 
force and diffusion coefficient of a matter wave bright soliton immersed in a thermal cloud. 
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Solitons are localized waves that propagate without 
spreading and attenuation. They appear from classical 
systems like ocean waves to optics and quantum systems 
like Bose-Einstein condensates (BEC) of atomic gases. 
A BEC of a dilute atomic gas with attractive two-body 
interactions in three dimensions (3D) is unstable and 
collapses -"- n one dimension (ID), however, a BEC 
with attractive interaction is stable against collapse and 
forms a self-bound particle-like object known as a bright 
soliton. Recently, bright solitons of Bose-condensed 7 Li 
atoms were observed in quasi- ID waveguide traps at Rice 
University and at ENS in Paris 0. 

One of the most important features of solitons is the 
dissipationless motion over long distances. Because of 
this property, optical solitons have important applica- 
tions in transatlantic fiber-optic communication systems 
. In this Letter we discuss how dissipative effects in the 
motion of a soliton in a thermal cloud can arise due to the 
3D nature of the BEC in a tight waveguide. Although 
the dynamics is strictly one-dimensional, the transverse 
extent of the mean field generates extra nonlinear terms 
in the effective ID equation, as first discussed in Ref. 
This formalism can also be applied to other systems de- 
scribed by a cubic nonlinear Schrodinger equation per- 
turbed by a quintic nonlinear term. An overview of re- 
lated work in nonlinear optics can be found in Ref. 0. 

In this Letter we consider the scattering of quasipar- 
ticles by the quasi- ID matter- wave soliton similar to 
linear waves scattering on breathers 0. Muryshev et 
al. considered the interaction of quasiparticles with dark 
BEC solitons and conjectured that these lead to accel- 
eration and eventually disintegration of the soliton in a 
thermal environment 0. Following a similar line of ar- 
guments, we show that the quasiparticles scattering on a 
bright soliton have a finite probability of reflection only 
due to the extra nonlinearity, which finally gives rise to 
dissipative effects. The bright soliton experiences friction 
and diffusive motion in a thermal cloud but maintains its 
integrity in contrast to dark solitons which disintegrate. 

A BEC in a waveguide with a harmonic transverse con- 



finement is well described by the Gross-Pitaevskii (GP) 
equation: 
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where ip is the macroscopic wavefunction of the conden- 
sate, lu is the frequency of transverse trapping poten- 
tial, and a is the 3D scattering length. In the quasi-lD 
limit the effective dynamics of the system takes place 
along the free axis (x-axis) without exciting the trans- 
verse modes. The quasi-lD limit can be achieved when 
the mean field interaction is smaller than the radial exci- 
tation frequency, 47r?i 2 |a||'0| 2 /rn < fuv. Aiming at an adi- 
abatic separation of slow longitudinal and fast transverse 
motion we can write the full 3D wave function assum- 
ing cylindrical symmetry as, ip(r,t) = (j)(x,t)x(p,x,t). 
Here, <p is the ID (longitudinal) wavefunction and x is 
the radial wavefunction with the normalisation conven- 
tion J \x\ 2 2Trpdp = 1 and J \cf>\ 2 dx — N, where N is 
the number of bosons in the system. In the adiabatic or 
Born-Oppenheimer approximation we now assume that 
the radial wavefunction \ depends only weakly on the 
slow variables x and t and their derivatives of x can be ne- 
glected. For this assumption to be correct, the time scale 
of dynamics should be longer than the inverse of the tran- 
verse frequency to and the relevant length scale should be 
significantly larger than the transverse oscillator length 
I = yjh/ (muj). After substituting the ansatz for the 
wavefunction into Eq. Q and neglecting the derivatives 
of x with respect to x and t, we obtain the following adi- 
abatically decoupled equations for the longitudinal and 
the transverse wavefunctions: 
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where we have introduced the transverse chemical po- 
tential p, which has to be found from the ground state 
solution of Eq. J2J as a function of the linear density 



2 



n(x,i) = |<^(a;,i)| 2 . A simple scaling argument shows 
that jl = hu>f(an), where /(■) is a dimensionless func- 
tion, which has been computed numerically in Ref. 0. 
Physical solutions of Eq. J3J) are found only if —an < 0.47 
H3- CiZ3l - otherwise transverse collapse occurs ^l|. I* 1 the 
following we will be interested in the quasi- ID regime of 
small an and expand p,(an) in a power series. 

In the quasi-lD limit, when \a\n <§; 0.47, the radial 
wavefunction \ wm be close to the ground state of the 
2D harmonic oscillator with a Gaussian profile. We 
can expand \ m terms of the radial eigenmodes tp u (p), 
X(p, x) = (fo{p)+J2 L/ C v {x)<Pu{p)- The coefficients C„ are 
small and can be calculated perturbatively. The trans- 
verse chemical potential jl can be obtained by using sec- 
ond order perturbation theory: 



turbation Hi are given by 



fi = huj + gn — gin + 



(4) 



where g = 2ah\o and gi = 241n(4/3)a 2 fiw. A correction 
to the ID coupling constant g beyond the GP approach 
presented here has been found in Ref. |l2j . The constant 
92 was calculated first in Ref. || and corrections beyond 
GP can be obtained by the self-consistent Hartree Fock 
Bogoliubov approach of Ref. 13]. We obtain the follow- 
ing effective equation describing the condensate in the 
quasi- ID limit: 



[-(h 2 /2m)d 2 x +g\^\ 2 -gM 4 ] ( j ) = fi ( j ) . 



(5) 



This is a nonlinear Schrodinger equation with a cubic 
and a quintic nonlinearity, as used before in Ref. [f| . The 
possibility of collapse is inherent in this equation as the 
quintic nonlinear term is attractive. An estimate from 
the 3D GP equation 0) gives stability of a single soliton 
solution if N\a\/l < 0.627 is fulfilled Q. 

Without the extra nonlinearity associated with g 2 , 
Eq. (J5J is integrable. For attractive interactions at 
a < 0, the bosons form a self-bound particle- like 
state known as a bright soliton with the wavefunction 
4>(x) = y N / 2b sech(x / b) and the chemical potential 
fx = -h 2 /2mb 2 : where b = l 2 /(N\a\). We notice that 
for a weak soliton parameter AT|a|/Z < 1, the system be- 
comes quasi-one-dimensional (b > I) . 

A soliton can be considered as a macroscopic particle 
of mass mJV, moving in the bath of thermal excitations. 
Dissipative motion of the soliton arises due to the scat- 
tering of thermal atoms. Here we consider the interaction 
of thermally excited particles with the soliton within the 
Bogoliubov formalism 0, 



[H + H 1 W = e(k)iP 
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where, ip — (u, v) is a two component vector of particle 
(it) and hole (v) amplitudes, and e is the quasiparticle 
energy. The unperturbed Hamiltonian Hq and the per- 
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where Vi = 2g\<j)\ 2 - 3. 92 |^| 4 and V 2 = g<j) 2 - 2g 2 \<j>\ 2 tj> 2 . 
The scattering states have energy e(fc) = ^m~ + In 
one dimension, neglecting the extra nonlinearity (32 = 0), 
we obtain the exact solution of the scattering states: 



Uk =A(k) [kb + i tanh(ar/&)] 
v k =A(k)scch 2 {x/b)e ikx , 



2 ikx 



(9) 
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where A{k) = l/(k 2 b 2 — 1) is a normalisation constant. 
The transmittance is given by 



t= (kb + if I '{kb - if 



(11) 



and the transmission probability is \t\ 2 = 1. Hence, the 
quasiparticles scatter without reflection on the soliton 
but only aquire a phase shift and a time delay in the 
scattering process. Reflectionless scattering on a soli- 
ton in the integrable nonlinear Schrodinger equation 
(with <7 2 = 0) is a well-known result of mathematical soli- 
ton theory and is also found in an exact solution of the 
quantum many-body model in the limit of large particle 
number [l^. In the quasi-lD limit, the soliton thus be- 
comes transparent and exhibits dissipationless motion in 
a thermal cloud. 

Now we consider the scattering problem of quasiparti- 
cles in the presence of an extra nonlinearity that breaks 
the integrability. Assuming that the coupling constant 
g2 of the extra nonlinear term is small, we can solve the 
scattering problem using Green's function techniques. In 
order to solve Eq. © for the particle amplitude u we 
construct the Green's function for the u component of 
H satisfying 

[-{h 2 /2m)dl — fi — e] Gx{x - x 1 ) = -5{x - x'), (12) 

which is given by G\(x — x') = (m/h 2 k) siri(k\x — x'\). 
Since the potential is symmetric, the scattering states 
can be constructed with even or odd symmetry. The 
Lippmann-Schwinger equation for the particle channel 
can be written as 



s /o = u ° e /o + J Gi(x - x l )V 1 (x l )u e/o (x')dx' 

+ [ G^x-x'^x'^ix^dx', (13) 



where u e / denotes even [odd] wave functions of the par- 
ticle states and u° e = cos(kx), [it,® = sin(fex)]. The most 
general wave function can be constructed from even and 
odd eigenstates: u k = Aul + BuZ- Asymptotically this 
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wave function becomes linx^ 
lim 



i^oo uj; = te lkx where \t\ 2 and R = \r\ 2 are the 
transmission and the reflection coefficient, respectively. 
We obtain R(k) by solving Eqs. EJ numerically and 
matching with the asymptotic solutions, see Fig. 
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FIG. 1: Reflection coefficient of a soliton as a function of 
momentum, for N\a\/l — 0.1, 0.2, 0.3, 0.4 (from smaller to 
larger values of reflection coefficient). The dotted line shows 
the analytical estimate 1141 for N\a\/l — 0.2. 

An analytical estimate of the reflection coefficient 
can be obtained from Eq. Ijl3|l by approximating <f> — 
\J N/2bsech(x/b) and u e / D and v e / a with the properly 
symmetrised solutions I10|) . This approximation be- 
comes exact for gi = and relies on gi being a small 
parameter. The reflection coefficient is given by R — \r\ 2 
and 

r(k) = -i -m ( 14 ) 



where the terms ij 



T2A(k) 



(J+ -<)(/_ +*)' 
and 2_ are given by 

J\r2l„l2 m 

fc6 + 61n(4/3)- 
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with Q±{x) = l/3 + :r 2 ±(l + x 2 ) 2 7rx/[3sinh(7rx)]. By us- 
ing the Lippmann-Schwinger formalism instead of a sim- 
ple Born approximation we obtain the correct limiting 
behaviour for small k where R — ► 1. Total reflection is ex- 
pected whenever the special resonant conditions leading 
to reflectionless scattering at 172 — are broken, as k — > 
implies a vanishing group velocity de/d(hk) — hk/m. 
This case is very different from phonons scattering on 
a perturbed dark soliton, which becomes transparent for 
small k as found in Ref. The approximation 1141) 
reproduces the qualitative features but slightly overesti- 
mates the exact values of R as seen in Fig. ^ 

The reflection coefficient R is a function of dimension- 
less momentum kb and the soliton parameter N\a\/l. 
In the dissipative dynamics of a macroscopic object 
like a soliton, the microscopic parameter N\a\/l enters 
through the reflection coefficient of the quasiparticlcs. 
Once we know the interaction of particles with a soli- 
ton from the microscopic theory, we can describe its mo- 
tion in the bath of thermal particles at a given tem- 
perature. A bright soliton is a mesoscopic object with 



mass mN, and its dynamics is governed by classical 
motion. Therefore, we can define a phase space distri- 
bution function of soliton's center of mass coordinate 
fiPiQit). When the soliton follows the classical trajec- 
tory then the distribution function takes a simple form 
f(p,<l,t) = S(p-p(t))5(q-q(t)), where p(t),q(t) are clas- 
sical phase space trajectories. In the presence of a bath 
of thermal atoms, the atoms impart a momentum to the 
soliton in the scattering process. While the soliton is at 
rest, the force imparted on the soliton cancels on the av- 
erage but, nevertheless, the stochastic nature of the force 
introduces a diffusive motion of the soliton. For a moving 
soliton, the average force imparted by the thermal par- 
ticles does not vanish and gives rise to a frictional force 
on the soliton. To include the dissipative effects in the 
soliton's motion we write down the kinetic equation for 
the phase space distribution function of the soliton |16( : 



dt dp \ dq J dq 



dH 
dp 



f =1 



coll ) 



(15) 



where, for small momentum transfer, the collision inte- 
gral Icon can be written as 



T ° 
-icon — Tr- 
ap 



(16) 



The terms A and B gives rise to friction and diffusion 
of the soliton respectively. The frictional force A can be 
computed from the following expression, 



A = 



dk 
2^ 



(-2hk)R(k) 



de(k) 



hdk 



N(E,k B T) (17) 



where N(E, ksT) describes the distribution of thermal 
particles in the frame of the moving soliton with veloc- 
ity v and the energy E takes the value E(k) — (hk — 
mv) 2 /2m. In each collision, the particle with momen- 
tum k has a probability R to reflect back and transfer 
the momentum —2hk to the soliton. This momentum 
transfer multiplied with the number of particles coming 
from each direction per unit time gives rise to a frictional 
force. When the soliton is at rest, the momentum transfer 
on each direction cancels on the average and as a result 
the friction vanishes. 

At finite temperatures the thermal atoms are dis- 
tributed according to the rules of quantum statistics. Al- 
though thermal equilibrium may be reached in an exter- 
nal trap j the subtle conditions of equilibrium are not 
necessarily fulfilled in a dynamical experimental situa- 
tion. Here, we consider the motion of a soliton relative 
to a significantly warmer thermal cloud of atoms. Wc 
thus can assume a classical Boltzmann distribution of 
thermal atoms, N(E, ksT) ~ exp(—E/kBT). Dissipa- 
tive effects of the soliton can be enhanced by increasing 
the density or the temperature of the thermal cloud. We 
consider the situation where 10 4 thermal particles are 
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confined within a length L = 50b (« 70/im for ENS soli- 
ton with b = 1.4^m), with a density of the thermal gas 
of nb = 200 (« 10 12 /cm 3 ) and the velocity distribution 
of the thermal particles being controlled by changing the 
temperature. Within a certain range of the soliton ve- 
locity the frictional force increases linearly with velocity 
as seen in Fig. [21 When the velocity is increased further, 
nonlinear effects take over and the force decreases. 




0.5 1 1.5 0.5 1 1.5 

k B T/||l| V/C s 

FIG. 2: a) Logarithm of the friction force A in units of fuui/l as 
a function of temperature for soliton velocity V = O.lh/mb, 
density of thermal gas nb = 200, and soliton parameter 
N\a\/l = 0.2 (solid line), 0.3 (dashed line), 0.4 (dotted line), 
b) Friction force as a function of velocity, for 7V|a,|/Z = 0.3, 
fcsT/|/x| = 0.5 (units and parameters as Fig. 2a). 

Now we can calculate the diffusion parameter of the 
transport equation: 



B = J ^2(hkfR(k)\^-\N(E,k B T) 



(18) 



This term describes the velocity fluctuations of the soli- 
ton and gives rise to a diffusion in the momentum space. 
A graph is shown in Fig. [21 
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FIG. 3: Logarithm of the diffusion coefficient B (in units of 
hmu! 2 ) of a soliton with zero average velocity as a function 
of temperature for the soliton parameters N\a\/l = 0.2 (solid 
line), 0.3 (dashed line), 0.4 (dot-dashed line) and a density of 
the thermal gas of nb = 200. 

So far we considered only the low-energy elastic scat- 
tering of thermal quasiparticles on the soliton. To restrict 
our discussion to the quasi-lD case, we neglected the 
higher energy radial excitations ~ ftio. As an additional 
effect, the soliton can radiate particles if the colliding 
quasiparticle has higher energy than the binding energy 



^| . Also nonlinear collective motion of the thermal cloud 
and the soliton is possible 0. However, the elastic scat- 
tering process discussed in this work will dominate if the 
condition kgT < < hoj is fulfilled. A tight radial 
trapping potential is suitable to avoid inelastic scatter- 
ing processes. In the ENS experiment 0, the oscillator 
length of radial confinement was I = 1.4/xm. For a soli- 
ton parameter N\a\/l — 0.4, the sound velocity at the 
center of the soliton becomes c s rs 2.5mm/s. If a soli- 
ton with N ~ 10 3 particles moves with a velocity 0.1c s , 
then it decelerates 5.01mm/s 2 due to the frictional force. 
Finally it stops after 0.05s, travelling a distance of 6/zm. 
The slowing down of a bright soliton can be observed 
experimentally by suitably manipulating the density and 
temperature of the thermal cloud. 

Due to the friction force, the momentum p of the soli- 
ton changes as dp/dt = —A(v). For small velocities, 
A = 7i> and the moving soliton stops after a time scale 
t = mN/j. Due to the diffusion process the energy of 
a resting soliton changes as E = (5/27) [1 — e -27 */ mJV ]. 
For t -C t, the energy of the soliton increases and finally 
it reaches a steady state with energy E = B/2-f. 

In conclusion, we have investigated the effects of a 
thermal environment on the dynamics of bright matter- 
wave solitons and have calculated the frictional force and 
diffusion coefficient in a microscopic approach. Friction 
and diffusion effects occur due to the deviation from the 
quasi-one-dimensional limit. Both of them are generally 
small and can be controlled by the parameters of the sys- 
tem if unattenuated propagation of solitons is desired. 
However, the parameters can be chosen such that the 
dissipative effects become accessible to experimental ob- 
servation with currently available techniques. 

We acknowledge enlightening discussions with 
G. Shlyapnikov who suggested this problem. 
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